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g Abstract 

■ A model for deeply virtual Compton scattering, based on analytical light-cone hadron wave func- 

tions is presented and studied at energies currently accessible at Jefferson Laboratory, DESY, and 
beyond. It is shown that perpendicular vector components play an important role at Q 2 < 10 GeV 2 
and that the meson-exchange diagrams are important at all energies. This could significantly im- 
pact the physical interpretation of the underlying hadronic amplitudes. 
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I. INTRODUCTION 



Our knowledge about the partonic content of baryons has for several decades been ob- 
tained from experiments in the deep inelastic scattering (DIS) region. At sufficiently high 
momentum transfer, the inclusive ep — > e'X reaction is very efficiently described in terms 
of parton distributions satisfying Bjorken scaling. The main reason for applicability of DIS 
to studies of hadron structure is that the process is factorizable into a hard part, calculable 
from perturbative quantum chromodynamics, and a soft part describing the quark-parton 
content of the hadron. This soft part gives the probability of finding a parton inside a proton 
with a certain fraction of the proton momentum. The soft part is usually parametrized in 
terms of parton distribution functions. 

In recent years, exclusive reactions like ep — > e'p'j (deeply virtual Compton scattering or 
DVCS) or ep — *> e'p'M (hard electroproduction of mesons), have been explored experimen- 
tally and theoretically, since they promise to provide further insight into the parton distribu- 
tions of hadrons. It turns out that even in certain exclusive reactions there is a factorization 
theorem that separates a hard, calculable part from a soft part Q . It is customary to define 
the soft interaction by introducing generalized parton distributions (GPD's) 

saw 

, that 

describe the transition amplitude for removing a parton with a certain momentum fraction 
and then putting it back with a different momentum fraction. The GPD's contain informa- 
tion about correlations between different parts of the proton wave function, or, equivalently, 
between partons at different locations inside the proton. This geometrical view has been 
explored by Burkardt, Diehl, Ralston and Pire, and others Jg. The calculations of DVCS 
have recently been extended to include twist-three effects |7(. The factorization of DVCS 
is conveniently represented by the handbag diagrams (Fig. [TJ, where the lower blobs are 
parametrized in terms of generalized parton distributions. The same diagrams have also 
been applied to wide angle Compton scattering and photon annihilation into hadron-anti- 
hadron pairs [sj]. 

In the Bjorken limit (Q 2 — > oo, Xb = Q 2 /^p-Q constant), there are four independent 
GPD's; H(x, (, t), E(x, C, i), H(x, (, t), and E(x,(,t), where x and ( are respectively the 
light-cone momentum fractions of the struck quark and real photon, while t = A 2 is the 
momentum transfer squared. Three physically different regions can be distinguished for x 
and (. The domain < ( < x < 1 (£ — 1 < x < 0) corresponds to the removal and return 



2 





p-k 



FIG. 1: The handbag diagrams, defining the kinematics and our notation. 



of a quark (antiquark) with momentum fractions x (£ — x) and x — ( (—%), respectively. 
In the remaining region < x < £, the photon scatters on a virtual quark-antiquark pair, 
extracted from the proton. The latter situation could be described as a meson exchange, but 

interpretation is to use light-cone perturbation theory j2j and regard the meson exchange 
as a higher Fock state of the proton splitting of a quark-anti-quark pair [10]. With this 
notation, ( — > xb (Bjorken x) in the limit Q 2 —>■ oo (A fixed). In the limit of forward 
scattering (DIS), the if's reduce to the quark density and quark helicity distributions, e.g. 
H(x, 0, 0) = q(x) and H(x, 0, 0) = Aq(x). The E's do not appear in DIS. They are unique 
to the off-forward exclusive processes and provide information not accessible through other 
means. Finally the x-integrated GPD's are related to the nucleon form factors |2]. 

Q 

An alternative approach to DVCS using light-cone quark wave functions [9J] was suggested 



in Refs. 



10 



111 ] . The authors have used this idea to calculate DVCS on an electron in QED 



for large Q 2 , assuming that the electron temporarily splits into a virtual electron- meson 
pair. 

There are ambitious efforts under way to measure DVCS (and the hard exclusive meson 

n n 

photo-production) at a number of facilities, in particular at DESY [12f and JLab [13j. Exper- 
iments using high-energy muon beams are planned for the COMPASS facility at CERN 14 1. 
At JLab energies, the competing Bremsstrahlung or Bethe-Heitler (BH) process is larger 
than DVCS. However, by carrying out interference measurements (e + /e~ beam charge asym- 
metry and various spin asymmetries), the BH amplitude cancels out and only a BHxDVCS 
interference remains. At DESY the energy is sufficiently large that DVCS becomes larger 
than BH, though the measurements have low statistics, insufficient to extract differential 



cross sections. However, both facilities are able to measure asymmetries for relatively low 
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beam energies and momentum transfers. In this paper we will investigate to what extent 
the 'leading' amplitude actually dominates in this kinematic region. In order to address this 
question, we employ a model using effective analytic quark-diquark wave functions. This 
model allows us to study the dynamics of the amplitudes excited at these low energies. Some 
preliminary results were reported in Ref. 16]. This model is similar to the one presented 
in Ref. jlOj , but we calculate DVCS on a proton and focus especially on the features at low 
Q 2 , and we keep higher- twist terms. 



II. FORMALISM 



In this paper, we treat the proton as a quark-diquark state. For simplicity, we drop 
the proton spin degrees of freedom. Furthermore we will only consider the valence quark 
flavors, i.e., we ignore any sea quark or higher Fock state contributions to the proton wave 
function. Thus, for a proton, the SU(3) operator A/v2 is replaced by the SU(2) (isospin) 
operator r/y/2. In future articles we will relax these constraints on our model, but this 
simplified model will suffice for the present purpose of studying the Q 2 dependence and 
analytic structure of the DVCS amplitudes. Our model proton state has the form 

. . f dxid 2 ki ± dx2d 2 k2± 16ii 3 6(x — X\ — X2)5 2 (ku_ + k 2 j_) 
\ P ' ~ J 16vr 3 16^ ^/xTx^ 

where Aj, fi, and q are the quark spin (helicity), flavor (isospin), and color indices. The 
proton state normalization is 

(p'\p) = 16n 3 6(x'-x)6 2 (p' ± -p ± ). (2) 

We start from the Fourier transform of the ^*p — > jpf amplitude (in light-cone coordi- 
nates) 

T ++ = -iJ d 4 ye^'-y(p'\TJ + (y)J+(0)\p), (3) 

where J + (y) = i , (y) r y + l5^4'(y) is the electromagnetic current, and p(p') and q(q') are 
respectively the four-momenta of the initial (final) hadron and photon. In light-cone time- 
ordered perturbation theory this expression can be expanded to give the five different one- 
loop diagrams (assuming light-cone gauge q + = 0) shown in Fig. |21 
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FIG. 2: The five different time-ordered 'handbag' diagrams. 

These diagrams are essentially the same as those used in Brodsky et al. ^(J for DVCS 
on an electron. The two approaches differ in that Brodsky et al. include the electron spin 
and treat the system as a fermion-vector composite, while we assume the diquark to be a 
fermion and approximate the proton as a scalar. In addition, we employ soft Gaussian form 
factors for our proton wave functions. The time-ordered diagrams of Fig. 121 are represented 
by the integrals 

{z, l±)<p(x, kj.) 
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In Eqs. EHHlwe define the kinematic variables z — (x — (,)/(! — Q, z' — (£ — x)/{l — x), and 



the relative momenta are defined as 



k_L + (1 - z)A ± 



1' 



■± 
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mi 
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k± + q± - -zq± 

k.^A, 
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m_L e 



-ki - -q ± 



(9) 



Here x = k + /p + > and £ = q' + /p + are the longitudinal momentum fractions of the struck 
quark and the real photon, the quantities m, rriR, and M are respectively the quark, remnant 
(diquark), and hadron masses, and A = q — q'. 

In the limit Q 2 — > oo, the denominators of Eqs. (@J) and (JJJ are proportional to a; — £ + ze 
and x + £ — it. These are the leading twist expressions of Ji 2]. In our model these 
leading-twist diagrams (a and d of Fig. EJ) are integrated over ( < x < 1 for which the 
Compton scattering occurs on a quark extracted from the proton. In principle, scattering 
on an anti-quark would have the same form, but it would require a different (higher Fock 
state) wave function. For the skew diagrams b, c, and e of Fig. the photon scatters on a 
quark- ant iquark pair split off from the proton. The proton wave function is represented by 
the analytic form 



agrees with the dipole form factor |17| for Q 2 < 1 GeV 2 . The skew diagrams require 
knowledge of the wave function <p(z', lj_) for the diquark splitting into a hadron and quark 
(the lower right-hand corner of diagrams b, c, and e of Fig. EJ) . The form of this wave function 
will eventually be restricted by exclusive data, but in this paper is arbitrarily chosen to be 
of the form of Eq. (|1U|). with tur and x replaced by M and z'. Since the relative phase of the 
two wave functions is not known, an arbitrary complex phase factor exp(«0) is introduced 
between the regular and skew diagrams. 




(10) 



where (3 = 0.6911 GeV is chosen such that 




(11) 
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The full T matrix could be related to the purely kinematic T ++ component by the Lorentz 

n 

structure T^ v oc p^n v + n^p v — ', derived by Ji |2|]. In this expression p and h are two 
lightlike four-vectors that project out the + and — components of Ji's formalism, and have 
been Lorentz transformed into our reference system. The Bethe-Heitler cross section is 
calculated from standard methods, using the dipole parametrization 17] of both the electric 
and magnetic proton form factors. Since we do not include the proton spin, the BHxDVCS 
interferences are not calculated. 



III. RESULTS 

The T ++ matrix elements have been calculated for present JLab and DESY kinematics, 
using standard Monte Carlo techniques. These calculations incorporate both the principal 
value (3?eT ++ ) and ^-function (9mT ++ ) parts. The angular distribution and Q 2 evolution 
of the real parts are plotted in Fig. El The laboratory angle # 7 y between the virtual and real 
photons is defined for in-plane kinematics such that it is positive for = and negative for 
= 180°, where is the azimuth angle between the final electron and the real photon, with 
q as the polar axis. Note that this angle convention is the reverse of the one in Ref. [16j |. 

The leading twist diagram (a of Fig. |2J is the largest for most scattering angles, with 
diagrams c and e giving contributions of roughly the same order. In our calculations, the 
contribution from the crossed diagram d is generally related to the leading-twist term by 
d ~ —a/10. This situation is quite different from DIS, where in general d <C a. This 
difference can be understood by realizing that DVCS is an exclusive process with an on- 
shell off-forward final photon (=>- SrnTd = 0, dteTd ^ 0), while DIS has an off-shell photon 
(SmTjj ^ 0), so that the imaginary part of the cut diagram is the only contribution to the 
DIS cross section. Since we are interested in the real part of diagram d (which is absent 
in DIS), and DIS uses the (negligible) imaginary part (which is absent in DVCS), the two 
reactions turn out to be quite different and do not measure the same physics. 

In the physical region of the DVCS process, diagrams a and b possess a cut on the real 
axis, while the other diagrams never get vanishing denominators for a massless final photon. 
Thus the imaginary part of the DVCS amplitude, which arises from the 5-function part 
of the propagators, has contributions from the first two diagrams (a and b) only. These 
imaginary parts are not plotted in Fig. |H] but are included in the calculation of the cross 
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FIG. 3: The KeT ++ matrix elements as functions of the laboratory angle between the virtual 
and the real photon for in-plane scattering. The labels correspond to the ones of Fig. [21 (the 
lower solid line always corresponds to diagram d). The calculations are for xb = 0.35 and Q 2 = 
2, 4, 10, 20 GeV 2 from top to bottom, using an arbitrarily assigned phase angle <p = 0°. 
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section. 

Because of its two hard propagators, diagram b is suppressed in the high Q 2 limit [lOj, and 
in our calculation this suppression does becomes significant, but only for large Q 2 > 10 GeV 2 , 
i. e., well above present JLab energies. On the other hand, diagrams c and e remain important 
at all energies and could be interpreted as contributing to meson exchange mechanisms, since 
the quark-antiquark pair could form an intermediate (off-shell) meson. The skew diagrams 
b, c, and e yield non-symmetric angular distributions, reflecting the shifted denominators in 
these cases. 

The difference in analytic structure also explains why diagram d gives a smaller contri- 
bution than diagram a, since the Q 2 terms add up in the denominator of d, while they have 
opposite signs for a. This feature is closely related to the behavior of the two leading-twist 
propagators of Ji [2] . These two terms have the form l/(x — £ + ie) and l/(x + £ — ie), where 
x and £ are momentum fractions related to +p' + ) instead of p + . In this notation there 
is a cancellation for x = £ (scattering on quark) or x = — £ (antiquark). The important 
difference between the two approaches is that Ji neglects four- vector components, e.g., A± 
that do not give large scalars in the Bjorken limit, while we explicitly keep them since we 
are interested in the behavior at low and moderate Q 2 . 

In Fig.|3]our predictions for the differential cross sections for DVCS and BH are shown for 
a range of beam energies, and for various choices of the relative phase between the regular 
and skew diagrams. These plots are comparable to those published in yl for the kinematics 

n n n 

corresponding to the JLab |13j |. DESY [121 ] . and CERN |1J| experiments. Only for very 
large Q 2 is it possible to separate DVCS from BH when only unpolarized cross sections are 
available. The dip in the DVCS cross section at ~ 2° for zero phase angle reflects the fact 
that the real parts add up to zero there. If such a structure would be observed (at large 
beam energies), it could help constrain the phase between the regular and skew diagrams. 

The effect of the skew diagrams on the cross section is seen Fig. El When they are removed 
(lower graph, solid line), the dip in the cross section disappears, but the cross section for 
backward angles also decreases by roughly an order of magnitude. However, the cross section 
is still asymmetric. If in the remaining leading twist diagrams a and d we set = in 
the energy denominators (lower graph, dashed line), we end up with a narrower, almost 
symmetric distribution. At small angles this shift increases the cross section by 30-40% and 
for backward angles it decreases it by a factor of about two. 
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FIG. 4: Differential cross section for DVCS (solid line) and BH (right-most dotted line) at xb = 0.3, 
Q 2 = 2 GeV 2 , and beam energies of the JLab, DESY, and CERN experiments. In the second graph 
the phase between diagrams a, d and b, c, e is 0° (solid line), 90° (dashed line), 180° (dotted line), 
and 270° (dot-dashed line). 
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FIG. 5: Effect of excluding the skew diagrams, for energy = 200 GeV. Upper graph: skew 
diagrams included. Lower graph, only the leading twist diagrams a and d are included. Lower 
graph, dashed line: only leading twist diagrams, and setting kj_ = in the denominators (dashed 
line). Note that the two plots use different linear scales. 
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The skewed behavior of the angular distribution is hence caused by a combination of 
two effects; the backward shifted meson exchange (skew) diagrams, and the presence of 
perpendicular components in the propagators. 

IV. CONCLUSIONS AND OUTLOOK 

We have introduced a model for DVCS, using simple analytic quark wave functions and 
retaining all components of every four- vector in the scattering process. This model enables us 
to investigate the five single-loop diagrams and their relative importance for DVCS at JLab 
and HERA energies. Our results indicate that at present JLab energies (Q 2 < 4 GeV 2 ), 
all of the diagrams of Fig. El (except possibly the crossed diagram d) contribute to the 
DVCS cross sections. In particular it is necessary to include diagram b, despite its two hard 
propagators, since the suppression of this amplitude is significant only for very large Q 2 . 
Even at kinematics appropriate for the upgraded (12 GeV) JLab, the maximal Q 2 ~ 6 GeV 2 
is too small to suppress diagram b. Other higher-order diagrams might need to be considered 
as well. At higher values Q 2 > 10 GeV 2 , the process is completely dominated by the handbag 
diagrams with one hard propagator. The crossed, w-channel diagram d of Fig. ^contributes 
on the 10% level. This should be contrasted with DIS, where this amplitude is very small 
and negligible. In our calculations the skew diagrams c and e are of the same order as the 
regular handbag a for all Q 2 . Consequently, meson exchange mechanisms are likely to be 
important for the understanding of DVCS. 

In the case of DIS, where cross sections are given by the imaginary part of the forward 
7*A — > 7*A amplitude, in the scaling region the t-channel meson exchange is relevant only 
at small xb- In the limit of large s and low t, where s = Q 2 (l — xb)/xb, the amplitude is 
expected to be dominated by the right-most singularity in the complex-/ plane. This leads 
to an amplitude proportional to P n (Q 2 )(Q /xbttin)""^ ■ Due to the large photon virtuality, 
however, the Regge residues /3 n (Q 2 ) are suppressed for high spin states (3 n {Q 2 ) ~ (1/Q) a ^°\ 
This simply follows from the t-channel sub-process amplitudes, 7*(<z)7*((/) — > n>i, where 
I = a(m 2 n ) 

M Q) ~ AW - n) <x /« • q) ( ,_ ( ^ )/2)2 «(!)'. (12) 

Thus, for DIS processes the Regge limit only applies at low i B , while at finite xb an 
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alternative (parton) description becomes more efficient. In the case of DVCS, however, 
the presence of the real photon in the exit channel invalidates the above argument, since the 
residue functions, (3 n {Q) ~ ^4(7*7 — > n{) oc 1/Q 2 are not suppressed for high spin states [18]. 
In other words, the presence of the soft q —> Nq fragmentation described by the bottom-right 
part of diagrams c and e masks the simple interpretation of the DVCS process in terms of 
partons originating from the nucleon alone. 

We have shown that the perpendicular components of four vectors significantly alters the 
angular distribution of DVCS, and makes it skew toward backward angles. This effect is 
further enhanced by the inclusion of the skew diagrams. 

This work will be extended to include a full treatment of the proton spin. This will 
allow us to calculate the interference of DVCS with the Bethe-Heitler process and hence to 
calculate spin asymmetries, which are needed for a thorough comparison with data. With 
this extended model, we will again test the validity of commonly applied approximations, 
and we will evaluate the effectiveness of various proton wave functions. We also intend to 
calculate meson photoproduction, with special consideration of the meson poles that appear 
in the skew diagrams. 

In conclusion, our calculations show that higher twist effects should be important in 
DVCS processes particularly at the relatively low energies currently accessible at JLab. In 
this kinematic region, higher twist amplitudes need to be carefully studied and calculated. 
In addition, a substantial part of the intuition we bring from DIS processes may no longer 
be valid in evaluating DVCS processes. 
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